Spin filtering by a periodic nanospintronic device 
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For a linear chain of diamond-like elements, we show that the Rashba spin-orbit interaction (which 
can be tuned by a perpendicular gate voltage) and the Aharonov-Bohm flux (due to a perpendicular 
magnetic field) can combine to select only one propagating ballistic mode, for which the electronic 
spins are fully polarized along a direction that can be controlled by the electric and magnetic fields 
and by the electron energy. All the other modes are evanescent. For a wide range of parameters, 
this chain can serve as a spin filter. 

PACS numbers: 71.70.Ej, 72.25.-b, 73.23.Ad 



I. INTRODUCTION 

In addition to their charge, electrons also carry a spin, 
which is the quantum relativistic source for the electron's 
intrinsic magnetic moment. Future device technology 
and quantum information processing may be based on 
spintronicSfi where one manipulates the electron's spin 
(and not only its charge). One major aim of spintronics 
is to build mesoscopic spin valves (or spin filters), which 
generate a tunable spin-polarized current out of unpo- 
larized sources. Much recent research aims to achieve 
this goal by using narrow-gap semiconductor heterostruc- 
tures, where the spins arc subject to the Rashba^ spin- 
orbit interaction (SOI): in a two-dimensional electron gas 
confined by an asymmetric potential well, the strength of 
this SOI can be varied by an electric field perpendicular 
to the plane in which the electrons move 4 An early pro- 
posal of a spin field-effect transistor— used the Rashba 
SOI to control the spin precession of electrons moving 
in quasi-one-dimensional wires. Placed between two fer- 
romagnets, the transport of polarized electrons through 
such a semiconductor could be regulated by the electric 
field. However, such devices are difficult to make, due to 
the metal-semiconductor conductivity mismatch. 

Some of the most striking quantum effects arise due 
to interference, which is best demonstrated in quan- 
tum networks containing loops. Indeed, interference 
due to the Rashba SOI has been measured on a 
nanolithographically-defined square loop array Here we 
discuss the possibility to construct a spin filter from 
such loops. Recently, several groups proposed spin fil- 
ters based on a single loop, subject to both an electric 
and a magnetic [Aharonov-Bohm (AB)^] perpendicular 
fieldsiiiSii^ However, such devices produce a full polariza- 
tion of the outgoing electrons only for special values of 
the two fields. In the present paper we consider a chain of 
such loops, as shown in Fig. [TJ The effects of the Rashba 
SOI on the spectrum of the diamond chain of Fig. [l]was 
studied by Bercioux et al.i^ They found a strong varia- 



tion of the averaged (over energies) conductance with the 
strength of the SOI, which they associated with localiza- 
tion of the electron due to interference between different 
paths in each diamond. Later, this groupii found similar 
effects due to both the SOI and an AB flux. However, 
the possibility to use such networks to achieve spin fil- 
tering has not been considered. As we show below, the 
polarization of the outgoing electrons depends on the en- 
ergy. Therefore, averaging over energies mixes different 
polarization directions and eliminates the possibility of 
obtaining full polarization. 
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FIG. 1: Chain of diamonds. 

We find that both the ballistic conductance and the 
spin polarization of the electrons going through the de- 
vice can be sharply varied by an electric field (deter- 
mining the SOI^), a magnetic field (determining the AB 
phases of the orbital electronic wave functions) and the 
electrons' energy (set by the chemical potential in the 
source). Varying these three parameters, we find large 
parameter ranges where all the energy eigenstates of the 
device except one become evanescent and decay exponen- 
tiall y, f orming the localized states discussed in Refs. [13 
and I 111 . However, the electrons in the remaining single 
mode propagate with fully polarized spins. Thus, elec- 
trons which enter with arbitrary spins exit fully polar- 
ized. Since this polarization can be tuned by the param- 
eters, our system is an ideal spin filter. 

Section II outlines the tight binding model which we 
use for solving the Schrodinger equation on the periodic 
chain of diamonds. Section HI presents results for the 
ballistic conductance and for the polarization of the elec- 
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trons in the regions where they are fuUy polarized. Fi- 
nally, Sec. IV contains a discussion of our results, includ- 
ing a comparison with the case of a single diamond and 
a discussion of the application of our results to a finite 
chain. 



II. TIGHT BINDING MODEL 

With SOI, we need to solve for the two-component 
spinor at each point on the network. Bercioux et ali^t^ 
treated each bond of the network as a continuous one- 
dimensional (ID) wire. Having expressed the solutions 
along each bond in terms of the spinors of the nodes at 
its two ends, they used the Neumann boundary condi- 
tions at the nodes to derive discrete equations for the 
spinors at these nodes. As we discuss elsewhere,— these 
boundary conditions are sufficient but not necessary for 
current conservation at the nodes. A more systematic 
way to treat such network replaces each continuous bond 
bond by a discrete sequence of sites, and then studies 
the tight binding model for the wave functions on these 
sites (and on the original nodes). As the number of these 
intermediate sites increases, one has more sites per unit 
cell, and therefore one ends up with more energy bands 
for the solutions which contain waves moving along the 
main axis of the network [i.e. along the (1,1,0) direction 
in Fig. [1] . Qualitatively, we find that all these bands are 
similar to each other, and also similar to those found for 
the continuous network used in Refs. [13 and [111. There- 
fore, we choose to report here only on the simplest case, 
with no intermediate sites within the bonds. Thus, we 
treat a simple tight-binding model, with sites {u} only 
on the corners of the diamonds. The latter model could 
also describe a network of quantum dots or anti-dots, lo- 
cated at these nodes The stationary spinors with 
energy e, obey the Schrodinger equations, 

ih{d/dt)^u = = -jJ2Uuv-^v, (1) 

V 

where the sum is over the nearest- neighbor nodes {v}, J 
is the (real) hopping matrix element (in the absence of 
fields) and 

U^,^Ul^eMK^uv +^'u?)] (2) 

is a unitary 2x2 matrix, representing the phase factors 
due to the AB flux and to the SOI, (j)^.^ and (j)^^ respec- 
tively. For our structure, all bonds are in the xy— plane, 
and both the uniform magnetic field H = Hz and the po- 
tential asymmetry which creates the SOI are along the 
z-axis. As can be seen from Fig. [TJ the n'th unit cell 
contains three sites, and Eq. (jT]) reduces to equations for 
the related spinors, ^'0(71), ^bin) and ^c{n)- Choosing 
the edges of the diamonds along the x— and y— axes (see 
Fig. [T]), so that site a„ is located at r„ ~ {n,n,0)L [L 
is the length of each edge) , the unitary hopping matrices 



within the n'th diamond are given hy^ 

C/a» = C^cK") = e-"'^/'e— ^ (3) 

where cr is the vector of Pauli matrices, a = k^pL 
{kgQ measures the strength of the 'microscopic' SOI, 
{h/m)kgQcr X p) and = 2'!tHL^ /^q represents the 
AB phase associated with a single square diamond (here, 
$0 = hc/e is the flux unit; h is Planck's constant, c is the 
speed of light and e is the electron charge). Note that 
the dependence of U^f^{n) and of Ug^^{n) on n results from 
our choice of gauge for the vector potential. The net flux 
through each diamond is equal to </>, independent of n. 

For e = one encounters dispersionless modes, for 
which "i/ain) = 0. Since these solutions have zero ve- 
locity, and therefore carry no current, we ignore them in 
the following discussion. We next eliminate the spinors 
^'{,(77-) and "^dn) from the equations, and end up with 
effective one-dimensional equations, 

4A*„(n) =Vt^',(n-l)+V5'„(n + l), (4) 

with 4A = (e/ J)2 - 4 and 

V ^ Uab{n)Uac{n + 1) + Uac{n)Uab{n + 1) 

Unlike the individual Uuv^s, the 'renormalized' hopping 
matrix V is not unitary. This lack of unitarity reflects 
interference between the two paths in a diamond, which 
may decrease the current along the chain. 

In the following we concentrate on propagating waves, 

vI/,(n)=Ae''^"Xa('z), (6) 

where L = L-\/2 is the lattice constant of the diamond 
system along its axis (1, 1,0), the (real) wave- vector q is 
in the range —Tr/2 < qL < tt/2 and Xa is a normalized 
spinor (which depends on q). For such solutions, Eq. 
Q implies that Xa must obey the eigenvalue equation 
T^Xail) = Axa('z), with the 2x2 hermitian matrix 

An e"*«-^Vt 4- e*«-^V. (7) 

We next write 

n = A + B-(T, (8) 

with 

A = cos{qL) cos((/>/2), 
B = — cs sin((7Z) cos((/)/2) 
X (l,-l,-cot(gZ)tan(^/2)s/c), (9) 

where c = cos a, s = sin a. It follows that the spinor 
Xa {q) must be an eigenvector of the spin component along 
n = B/|B|: n • (TXa{q) = t^Xa{q), = ±1- Thus, 
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A ~ A + /i,|B|. Given A, this equation can be writ- 
ten as a quadratic equation in x ~ cos(gL). Denot- 
ing the solutions by X12, we end up with four solutions 
(7^2 = i arccos 2 ■ These solutions are propagating 
(evanescent) if q is real (complex). For each q one then 
has /i = (A — ^)/|B|, so that /.j is invariant under flipping 
the sign of q. 

Since Xail) is an eigenvector of n • cr, each solution 
with a given q is associated with a full polarization along 
the direction n, 

S= (Xa(9)|T|xa(g)) -Mn. (10) 

As usual for Rashba SOI, n is always perpendicular to the 
direction of motion along the axis of the diamond chain, 
(1,1,0). In the absence of an AB flux (i.e. (p = 0) n 
remains in the direction (1,-1,0). However, the orbital 
AB flux causes a rotation of the polarization axis towards 
the z— direction. Below we present results for Sz and for 
Sxy = {Sx — Sy)/\/2. Since n^, ^ (n^) is odd (even) in q, 
flipping the sign of q flips the sign of Sxy but not that of 
Sz. 

The probability current from site u to site v is 

I{u ^v)^ {2J/h)'^{^,,\Uuv\^.). (11) 

The current from site a„ to site a„+i on the diamond 
chain, equal to the sum of the currents from a„ to &„ 
and to c„, is thus found to be 

I{n^n + l) = -{2J^/he)^[{^a{n)\V\^a{n + 1))]. 

(12) 

For a single propagating solution of the form ([6]), Eqs. 
© and Uni) yield 

I{n -{U'^ /he)\A\^[sai{qL){cos{(t>/2)c^ 

+ liriz sin((/)/2)s^) 

-I- y,{nx — riy) cos{qL) cos(0/2)sc] . (13) 

It is easy to see that / flips sign with q. When we have 
only a pair of propagating modes, we thus concentrate 
on the one with / > 0. 



III. RESULTS 

Figure [2] shows the spectrum t{q) of the propagating 
solutions (real g's), for several values of and a. The left 
column shows results for = 0, similar to Ref. in- 
creasing a splits the energy band vertically, and changes 
its width. Thus, the SOI can turn propagating waves 
into evanescent ones, with complex q (our figures show 
only the solutions with real q). However, whenever the 
energy e allows for real values of q, there exist four such 
values, forming pairs which move in opposite directions 
and have opposite spins along (1,-1,0). 

The situation becomes more interesting when we have 
both the SOI and the AB flux. Adding only the latter 



(upper plot on the right hand side of Fig. ^ creates a gap 
(i.e. evanescent states) around e = 0. The degeneracy of 
the propagating solutions is not lifted, since the two spin 
directions have exactly the same energies. As seen in the 
right column in Fig. [21 increasing a at fixed </> = QAtt 
causes the splitting of each sub-band horizontally. 

We next discuss the ballistic conductance of our device, 
G. For an ideal conductor, this conductance is given by 
G ~ {e'^/h)g, where g is the number of right-moving (or 
left-moving) propagating modes at a given energyj ^^'^^i^^ 
This formula clearly applies for the infinite periodic chain 
of diamonds discussed here. Below we argue that the 
filtering effect which we find also survives for a finite 
chain, under certain conditions. As Fig. [5] shows, at 
a given energy e one can encounter zero, two or four 
propagating solutions. The number g can be read directly 
from Fig. [2| on the left hand side of this figure, the 
number of real g's (both left-moving and right-moving) 
is always zero or four, and thus ,9 = or 2. In contrast, 
the right hand side of Fig. [5] shows 0, 2 or 4 real g's, i.e. 
(7 = 0, 1 or 2, depending on the parameters e, (j) and a. 




FIG. 2: (Color online) The spectrum {q versus e) of the prop- 
agating solutions. Here, J = 1 and the wave vector q is in 
units of ttL. Left: (p = 0. Right: <j) = OAn. Top to bottom: 
a — 0, 0.27r, OAtt. The vertical lines indicate boundaries at 
which the number of propagating solutions changes. 

We next consider electrons coming with arbitrary spin 
directions from a reservoir at —00, with energy e equal 
to their chemical potential in that reservoir. For each 
electron, its spinor will become a combination of the 
eigenmodes of the problem inside the system. In fact, 
the same will happen to electrons which enter into a fi- 
nite but long chain from the left hand side: their spinor 
within the chain will become a similar combination of the 
four eigensolutions there, multiplied by some transmis- 
sion coefficients. When all four g's have non-zero imag- 
inary parts, all of these modes are evanescent, and the 
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wave function will decay to zero, resulting with zero cur- 
rent. In that case there are no propagating modes, and 
5 = 0. When all four q's are real, i.e. g = 2, the incom- 
ing wave function is a combination of two right-moving 
modes, and it has no definite spin. However, for g ~ 1 
the wave function of the right moving electron is a lin- 
ear combination of one propagating and one evanescent 
modes. The latter will decay, and the spinor will con- 
verge to that of the single propagating solution, which 
has a uniquely polarized spin, see Eq. (jlOp . Without the 
AB flux, we always had g = 2 or g = 0. For 7^ 0, 
we find regions of energy where g = 1. Figure [3] shows 
contour plots of g in the (p — a plane, for several values of 
e. As one can see, for energies e = — 1.2J and e = — 2.4J 
there are large regions where g = \. In these regions, 
the electron will have a well defined polarization, which 
depends only on e, 4> and a. 

We next consider specific cuts through these contour 
plots. Figure m shows results as a function of a for fixed 
energy e/J = —2.4 and AB phase </> = 0.297r. The plots 
show only the (17 = 1 or g = 2) right-moving modes 
(/ > 0). The other propagating modes have opposite 
signs for q, I and Sxy- The dashed curves represent the 
second mode, which arises only when g ~ 2. For our 
purposes, we concentrate on the regions where g = I, 
where one has only the dotted lines. The top plots show 
the solutions for q and the corresponding currents /. The 
bottom plots show the spin components Sxy and Sz ■ The 
variation of Sxy with a is striking: the spins of the prop- 
agating electrons switch the sign of their in-plane compo- 
nent with a small change of a near a = ±0.57r. Note also 
the flipping of Sxy as a crosses ±7r. This flipping per- 
sists as (j) increases, and the range with g = 2 near these 
points narrows. Figure [5] shows results as a function of (j), 
for the same energy, but at flxed SO strength a = O.TStt. 
Clearly, even a relatively small AB flux already yields a 
single right-moving propagating mode {g = 1) and there- 
fore fully polarized spins. At small (j), the polarization 
starts close to the (1,-1,0) direction, but it then rotates 
towards the z— direction as </> increases towards ±7r, and 
flips sign after crossing these points. 



FIG. 3: (Supplied separately) Contour plots of the ballistic 
conductance (in units of e^/h) in the 4> — a plane (the AB 
phase (j) and the SO strength a are in units of tt). The values 
<; = 0, 1, 2 are represented by dark, medium and bright areas. 
The number above each plot is the energy e (in units of J). 



IV. DISCUSSION 

Given the above analysis, we may compare our sys- 
tem with that of the single diamond, Ref. As we re- 
port elsewherefi^ the single diamond generates fully po- 
larized electrons, along a controllable direction, whenever 
cos^(0/2) and for any e. Although this condi- 
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FIG. 4: (Color online) Wave vectors q (in units of nL), cur- 
rents I, and spin components Sxy and Sz, for right-moving 
modes, as functions of the SO strength a (in units of tt), for 
e/J = —2.4 and (j) = 0.297r. For values of a at which g = 1, 
the figures show only one mode (dotted line). When g = 2, 
the figures show two modes (dotted and dashed lines). 
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FIG. 5: (Color online) Same as Fig. H] for e/J = -2.4 and 
for fixed SO strength a = 0.75%, as functions of 4> (in units 
of tt). 



tion is less restrictive than that given in Ref. |8|, it is 
still much more restrictive than the conditions we found 
above. The literature contains many other proposals for 
spin fllters, also based on the Rashba SOI. Usually, these 
give only a partial polarization. Some of these devices 
also require a large Zeeman field. In contrast, our filter 
can work at a relatively low (and fixed) magnetic field (as 
apparently desired technologically), so that the Zeeman 
energy is negligible. Note also that both a and <j> depend 
on the diamond size L, and therefore one can choose a ge- 
ometry which corresponds to the available ranges of the 
magnetic field and the microscopic Rashba parameters. 

In real experiments it is not realistic to use an infinite 
chain of diamonds. We now argue that under appro- 
priate conditions it is sufficient to use a finite chain, as 
long as it is longer than the decay lengths of the evanes- 
cent modes. For the electrons coming in from the left we 
don't need to worry about the details of the connection 
between the incoming lead and the chain: even if some 
of the electrons are reflected back into that lead, those 



5 



which arc transmitted into the chain will split into a sum 
of the four modes there, and when g = 1 wc still remain 
with fully polarized electrons (although their overall am- 
plitude may involve a transmission factor with magni- 
tude smaller than 1). The situation on the right hand 
end of the chain is more delicate. Here we should avoid 
reflections, since they may modify the outgoing spinors 
and change their polarization. A standard way to avoid 
reflections is to use adiabatic contacts. This is usually 
done for retaining the ballistic conductance of mesoscopic 
devices.— One way to avoid reflections is to have a large 
leakage to the ground near the exit channel, so that only 
a small fraction of electrons enter into the exit lead. 

For our filter to be useful, one also needs to measure 
the outgoing spins, or to relate the outgoing spin polar- 
ization to some measurement of a voltage or a current. 
This issue is common to many proposed filters, and it 
requires separate research. For the present purposes, we 
mention just a few possibilities. First, one can follow 
the original proposal of Datta and Das;^ and connect 
the right hand end of the device adiabatically to a ferro- 
magnetic lead, whose magnetization can be tuned. The 
outgoing current will decrease with the angle between the 



electron polarization and this magnetization. Second, to 
avoid connections to fcrromagncts, one can also connect 
our filter adiabatically to another such filter, with differ- 
ent parameters which may block the polarized electrons 
coming from the first filter. 

Another way to test the spin polarization, is to couple 
one of the a-nodes (Fig. [T]) to a side quantum dot, that is 
in a Pauli spin blockade region.— After a while, the side 
dot will capture one of the polarized electrons, and this 
will block the current (which contains electrons with the 
same polarization). Changing the parameters will then 
change the spin direction of the propagating electrons, 
and allow some current until the next blocking occurs. 

In conclusion, we propose a simple spin filter, which 
yields a full polarization over a broad range of param- 
eters. For given energy e and magnetic flux (j) (which 
need not be very large), the polarization of the outgoing 
electrons can be tuned by varying the electric field which 
determines the SOI strength a. 

We acknowledge discussions with Joe Imry. AA and 
OEW acknowledge the hospitality of NTT and of the 
ISSP, where this project started, and support from the 
ISF and from the DIP. 
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